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Abstract 

We address existence of moving gap solitons (traveling localized solutions) in the Gross-Pitaevskii 
equation with a small periodic potential. Moving gap solitons are approximated by the explicit 
localized solutions of the coupled-mode system. We show however that exponentially decaying 
traveling solutions of the Gross-Pitaevskii equation do not generally exist in the presence of a 
periodic potential due to bounded oscillatory tails ahead and behind the moving solitary waves. 
The oscillatory tails are not accounted in the coupled-mode formalism and are estimated by using 
techniques of spatial dynamics and local center-stable manifold reductions. Existence of bounded 
traveling solutions of the Gross-Pitaevskii equation with a single bump surrounded by oscillatory 
tails on a finite large interval of the spatial scale is proven by using these technique. We also show 
generality of oscillatory tails in other nonlinear equations with a periodic potential. 



1 Introduction 

Moving gap solitons are thought to be steadily traveling localized solutions of nonlinear partial differ- 
ential equations with spatially periodic coefficients. The name of gap solitons comes from the fact that 
parameters of stationary localized solutions reside in the spectral gap of the associated Schrodinger op- 
erator. Existence of stationary solutions can be studied by separation of the time and space variables 
and reduction of the problem to an elliptic semi-linear equation. Since the variables are not separable 
for traveling solutions, a little is known about existence of moving gap solitons. 

We address stationary and traveling localized solutions in the context of the Gross-Pitaevskii equation 
with an external periodic potential 

iEt = -E xx + V(x)E + a\E\ 2 E, (1.1) 

where E(x,t) : R x R + k C, V(x) : R i-» R, and a = ±1. The Gross-Pitaevskii equation (fT7TT> is 
derived for the mean-field amplitude of the Bose-Einstein condensate placed in the optical lattice trap 
V(x), where a is the normalized scattering length [7j. Stationary solutions of the Gross-Pitaevskii 
equation are found from the solutions of the nonlinear differential equation 

4>"(x) + U3<f>{x) = V{x)(f){x) + a\(j)(x)\ 2 (j)(x), (1.2) 
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where (ft(x) : R i— > C, uj G R, and the exact reduction E(x,t) = (fi(x)e lujt is used. Localized stationary 
solutions of the ODE problem (|1.2|) were proved to exist in [9] . 

Theorem 1 (Pankov, 05) Let V(x) be a real-valued, bounded, and periodic potential. Let uj be in 
a finite gap of the purely continuous spectrum of L = —d 2 + V(x) in L 2 (R). There exists a non- 
trivial weak solution <j){x) in ff 1 (R), which is (i) real-valued, (ii) continuous oni£R and (Hi) decays 
exponentially as \x\ — > oo. 

In |10j . we have obtained a more precise information on properties of the stationary solution (j){x) 
by working with a small potential V(x), when the spectrum of L exhibits a sequence of narrow gaps 
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bifurcating near resonance points uj = ui n = \, n G N. We have justified the use of the stationary 
coupled-mode equations which have been used in the physics literature [13] for explicit approximations 
of stationary gap solitons. 

In this paper, we shall investigate whether the time-dependent coupled-mode equations can be used 
for approximation of moving gap solitons in the framework of the Gross-Pitaevskii equation (jl.ip . The 
coupled-mode equations are derived formally when the potential is represented by V = eW(x), where 
e is small parameter and W(#) is a normalized e-independent function described as follows. 



Assumption 1 Let W(x) be a smooth 2n-periodic function with zero mean and symmetry W(— x) = 
W(x) on x G R. The Fourier series representation ofW{x) is 

W(x) = ^ W2 m e imx , suchthat ^ (1 + m 2 ) s \w 2m \ 2 < oo, Vs > ^ (1.3) 

m£Z m£Z 

where wq = and W2 m = w_2 m = W2 m , Vm G N. 



An asymptotic solution of the Gross-Pitaevskii equation (|1.1|) is represented in the form [U [T2] : 



E(x,t) 



,1/2 



a(ex,et)e 2 + 6(ex, et)e 



1 a .f 
2 



tl t _ , a /o, 

— +0(e 3 / 2 ), 



where the vector function (a, b) : R x 



i(aT + na x ) 
i(br — nbx) 



satisfies the coupled-mode system 

w 2n b + o-(\a\ 2 + 2\b\ 2 )a, 
w 2n a + a{2\a\ 2 + \b\ 2 )b, 



in slow variables X = ex and T = et. The system fll .5 j) admits a separation of variables [13j : 



a 



n + c 



1/4 



A(0e 



-i/inr 



n — c 



1/4 



i/mr 



a — r j \n + c / 

where |c| < n, the new independent coordinates (£, r) are given by the Lorentz transformation 

X -cT T-cX 



72' 



n 



and the new functions A(£) and -B(^) satisfy the coupled-mode system 



(n — c) (iA 1 — w 2n B) + fin(l — cn)A 
-(n + c) (iB' + w 2n A) + fin(l + cn)B 



{n + c)\A\ 2 + 2{n - c)\B\ 2 
2{n + c)\A\ 2 + (n - c)\B\ 2 



A, 
B. 



(1.4) 



(1.5) 



(1.6) 



(1.7) 



(1.8) 
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Since \ A\ 2 — \B\ 2 is constant in £ G M and the constant is zero for localized solutions at infinity |£| — > oo, 
we can further represent the localized solution in the form 

A = <f>(£)e i rtt), B = $(£)e*e®, (1.9) 

where the functions </> : R i— > C and ^:Kh>R are solutions of the first-order equations 

, nc(u(l-n 2 )-2a\4>\ 2 ) 
V = J^-^) > (l w) 

■Ml I M« 2 (l-c 2 ) , . (3n 2 -c 2 ), ,,2 , 1 } 

The second equation of the system (ll.lOp is closed on <^>(£) and the explicit localized solution for c / 
can be found from the corresponding solution for c = [TO]. For instance, if a = —1 and win > 0, the 
function 0(£) is found in the explicit form 



] 1{n 2 -c 2 ) 04^1 (Ln) 

^ 3n ~ C ) V^2n + MO COsh (V«4i _ MoC) ~ V™2n ~ MO SUlh ( \/ w 2n ~ Mo? 



n 2 (l— c 2 ) 

where uo = M a\ an d Mo < ^2™- In the case c = 0, we have /Lto = A* an d the condition kt < iU2n 
indicates that the frequency parameter uj = u> n + e/i of the stationary gap soliton with c = is chosen 
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inside the newly formed gap of the continuous spectrum near the bifurcation point uj n = 2 |-, n G N 
[10], such that 

Wn - ei«2n < W < + e^2n- (l-l 2 ) 

The exact solution (jl.lip can be extended easily to values a = — 1 and u>2n < 0. Given a localized 
solution for </>(£), we can integrate the first equation of the system (jl.lOp and obtain a linearly growing 
solution for </?(£): 

V = (n2 n _ C c2) (^(1 " n 2 )t - 2a j* m')\ 2 dA . (1.13) 

The trivial parameters of translations of solutions in £ and <p are set to zero in the explicit solutions 
(jl.llj) and (|1.13j) . such that the functions A(£) and -B(£) given by the parametrization (II. 9h satisfy the 
constraints A(£) = A(-g) and B(£) = £(-£)• 



Definition 1 JTie traveling solution of the coupled-mode system jil.5\) is said to be a reversible homo- 
clinic orbit if it decays to zero at infinity and satisfies the constraints A(£) = A{— £) and 5(£) = B(—£) 
in parametrization U.6\) - [T^\) . 



We study persistence of the traveling solution (jl.6() — (11. 13|) with c 7^ of the coupled-mode system 
(jl.5|) in the Gross-Pitaevskii equation (jl.ip . We show that the moving gap solitons have bounded 
oscillatory tails in the far-field profile of the scale e N+1 , which are small in amplitude of the order 
e N+i/2 £ or an y > l. These small oscillatory tails are not accounted in the coupled-mode system 
(|1.5p . The main result is formulated below. 

1—1 2 

Theorem 2 Let Assumption^ be satisfied. Fix n G N, such that wm 7^ 0. Let u> = ^- + eSl, suc/i £/ia£ 

|f2| < r2o = |u) 2 n|- y -^-^-- Let < c < n, such that n2 + c2 ^ Z', where U is a set of odd (even) numbers 
for odd (even) n. Fix N G N. For sufficiently small e, there are e-independent constants L > and 
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C > 0, such that there exists an infinite- dimensional, continuous family of traveling solutions of the 
Gross-Pitaevskii equation M-l\) in the form E(x,t) = e~ tujt ^(x,y), where y = x — ct and the function 
if}(x,y) is periodic (anti-periodic) function of x for even (odd) n, satisfying the reversibility constraint 
ip{x,y) = tp(x,—y), and the bound 



ip(x,y) 



,1/2 



a e (ey)e^ + b e (ey)e-^' < C Q e N+1 '\ Vx G R, Vy G [-L/e N+1 ,L/e 



iV+li 



(1.14) 



Here a e (Y) = a(Y) + 0(e) and b e (Y) = b(Y) + 0(e) on Y = ey G R are exponentially decaying solutions 
as \Y\ — ► oo, where a(Y) and b(Y) are solutions of the coupled-mode system il.5\) with Y = X — cT . 



Remark 1 

(a) The solution ip(x,y) is a bounded non-decaying function on a large finite interval 

ye[-L/e N+1 ,L/e N+1 ]cR 



but we do not claim that the solution ip(x, y) can be extended to a global bounded function on y G R. 

(b) Since the homoclinic orbit (a, b) of the coupled-mode system (II. 5h is single-humped, the traveling 
solution tp(x, y) is represented by a single bump surrounded by bounded oscillatory tails. 

(c) The solution (a e , b f ) is defined up to the terms of O(e^) and it satisfies an extended coupled-mode 
system with a unique reversible single- humped homoclinic orbit. 



Our work can be compared with three groups of papers. The first group covers rigorous justification 
of the validity of the coupled-mode system (jl.5p for the system of cubic Maxwell equations [I] and for 
the Klein-Fock equation with quadratic nonlinearity [12] . The bound on the error terms was proved 
for a finite time interval, which depends on e. By using this bound, one can see that the solution of 
the Gross-Pitaevskii equation (jl.ip behaves as a moving gap soliton of the coupled-mode system (II. 5ft 
during the initial time evolution in i7 1 (R) [4J or in C?(R) [12]. However, the error is not controlled on 
the entire time interval t G R since other effects such as radiation due to interactions of the moving 
gap soliton with the stationary periodic potential can destroy steady propagation of gap solitons. 

The second group of articles covers analysis of persistence of small-amplitude localized modulated 
pulses in nonlinear dispersive systems such as the Maxwell equations with periodic coefficients pQ, the 
nonlinear wave equation [5], and the quasilinear wave equation [6]. Methods of spatial dynamics were 
applied in these works to show that a local center manifold spanned by oscillatory modes destroys 
exponential localization of the modulated pulses along the directions of the slow stable and unstable 
manifolds. As a result, the modulating pulse solutions decay in the spatial dynamics to small-amplitude 
oscillatory disturbances in the far-field regions. 

The third group of papers addresses propagation of a moving solitary wave in a periodic potential 
V(x) of a large period (see review in [H]). An effective particle equation is derived from the focusing 
Gross-Pitaevskii equation (II. ip with a = — 1 by an heuristic asymptotic expansion. The particle 
equation describes a steady propagation of the moving solitary wave with u < 0, which corresponds 
to the semi-infinite gap of the periodic potential. Radiation effects appear beyond all orders of the 
asymptotic expansion. They have been incorporated in the asymptotic formalism by using perturbation 
theory based on the inverse scattering transform [8j. The same methods were applied to the finite- 
period and small-period potentials [11] . Unfortunately, this group of article does not connect individual 
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results in a complete rigorous theory of the time evolution of a solitary wave in a periodic potential, 
although it does gives a good intuition on what to expect from the time evolution. 

Our article is structured as follows. Section 2 reformulates the existence problem for moving gap 
solitons as the spatial dynamical system. Section 3 presents the Hamiltonian structure for the spatial 
dynamical system and normal coordinates of the Hamiltonian system. Section 4 describes a transfor- 
mation of the Hamiltonian system to the normal form and gives a proof of persistence of a reversible 
homoclinic orbit in the extended coupled-mode system. Section 5 presents a construction of a local 
center-saddle manifold which concludes the proof of Theorem [2j Section 6 discusses other models for 
moving gap solitons with oscillatory tails. 



2 Spatial dynamics formulation 

We look for traveling solutions of the Gross-Pitaevskii equation (jl.ip in the form 

E(x, t) = e~ tujt i;(x, y), y = x - ct, (2.1) 

where a; is a parameter of gap solitons and the coordinates (x, y) are linearly independent if c 7^ 0. 
For simplicity, we only consider the case c > 0. The envelope function ip(x,y) satisfies the partial 
differential equation 

(w - icd y + d 2 x + 2d x d y + fig) i>{x, y) = eW(x)4>(x, y) + a\4>{x, y)\ 2 4>{x, y). (2.2) 

At this stage, the equation (|2.2I) is equivalent to the original equation if c 7^ 0. We shall however 
specify the class of functions t/j(x, y) to accommodate the moving gap solitons according to their leading- 
order representation given by ()1.4|) . (jl.6p . and (jl.7p . In particular, we consider either periodic (for even 
n) or anti-periodic (for odd n) functions ip(x, y) in variable x E [0, 2ir] and look for decaying or bounded 
solutions ip(x,y) in variable i/GR. Such solutions can be described by using the formalism of spatial 
dynamical systems [5j [6]. We make use of the periodic or anti-periodic conditions in variable x and 
represent the solution ip(x,y) in the form 

^(x, y) = E ^rn(y)e> x , Vm(y) = n~~/ = r V)e~ imx dx, m E Z', (2.3) 

where the factor yfi is used for the convenience and the set Z' contains even numbers if ^(x, y) is periodic 
in x and odd numbers if 1(1 (x, y) is anti-periodic in x. The series representation (j2.3|) transforms the 
PDE system (|2.2p to the nonlinear system of coupled ODEs 



mieZ' 



ra-mi — m2 

(y), m E Z'. (2.4) 



The left-hand-side of the system (|2.4p represents a linearized system at the zero solution for e = 0. 
Since the linearized system at e = has a diagonal structure on m E Z', its solutions are given by the 
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eigenmodes Vw(?/) = eKmV ^m,m' with m, m' G Z', where the values of K m are determined by the roots 
of quadratic equations 

2 

K = K m : k + i[m — c) k + u — = 0, m£S, (2-5) 

2 

The zero root k = exists if and only if uj = w n = ^- for any fixed n G Z. The zero root has multiplicity 
two for m = ±n if n G N and c 7^ ±n. The special value uj = ui n corresponds to the bifurcation of 
periodic or anti-periodic solutions as well as of the stationary gap solitons with c = [10]. We note that 
the values of n determine the choice for the set Z': it includes even (odd) numbers if n is even (odd). 
We shall hence focus on the bifurcation case uj = u n , when the two roots of the quadratic equations 
(I2.5h are represented explicitly as follows 



± i(c — m) ± V 2cm — n 2 — c 2 , 

r o . o " 

When m > uiq 



uj = uj n : k = K m = , m G Z . (2-6) 



n z +c z 



2c 



, where the notation [a]' denotes the integer part of the number a in the set 
Z', all roots are complex-valued with Re(ft^) = ±^y/2cm — n 2 — c 2 7^ and Im(/^j) = ^y^. When 
w < mo, al roots k are purely imaginary with = z/j^ and = ( c ~ m ) =l=v/ ^ 2 + c2 ~ 2c ^ j 

Lemma 1 Let uj = u) n , n G N and c > 0, such that n ^ c ^ Z'. Then, 

(i) The phase space of the linearized system \2.J$ at the zero solution for e = decomposes into a 
direct sum of subspaces E s © E u © E c+ © E c , where 

The subspace E^ consists of the eigenspace associated with the m-th Fourier component of the 
solution \2. 31) corresponding to the root (KM). 

(ii) The zero root k = is semi-simple of multiplicity two. The purely imaginary roots k G iM. are 
semi-simple of the maximal multiplicity three. All other roots k G C are simple. 

Proof. It follows from the quadratic equation (|2,5p that a root k = kq is double if kq = ^ c 2 , which 
implies that m = - ^ ■ Under the non-degeneracy constraint n ^ £ an roots k are semi-simple. 

n 2 +c 2 



When m > tuq 



, all roots are complex-valued and simple, such that E s and E u are stable 



2c 

and unstable manifolds of the linearized system at the zero solution for e = 0. 

When m < mo, all roots k are purely imaginary, such that E c+ © E c is a center manifold of the 
linearized system for e = 0. It is obvious that k+ increases as m decreases, while k~ decreases for 

r 2 i' 

m\ < m < mo and increases for m < mi as m decreases, where m\ = y- . Therefore, the purely 
imaginary roots may have the maximal multiplicity three. The zero eigenvalue has however multiplicity 
two since the two modes m = n and m = —n have simple zero eigenvalues and other modes have no 
zero eigenvalues. □ 

Lemma 2 Let uj = u n> n G N and < c < n. If c is irrational, all non-zero roots k of the quadratic 
equations $2.5\) are simple. 
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Proof. By Lemma [TJ only imaginary roots k can be semi-simple. Let two roots k coincide for m < ttiq 
and I < niQ. Then, k = —i^^- and (m, I) satisfies the equation 

(m - I) 2 + 4c(m + I) - 4n 2 = 0, m<m ,l<m . (2.8) 

If c is irrational, equation (|2.8p has no solutions for integers m and /. Therefore, all non-zero roots k 
are simple. □ 

Example 1 Figure 1(a) illustrates the distribution of imaginary roots = ik^, fn < mo for c = 
and n = 1. Although all imaginary roots are simple for this (irrational) value of c, the purely imaginary 
roots k can approach to each other arbitrarily close. Figure 1(b) shows a similar distribution for c=\ 
and n = 1. It follows from equation (|2.8p that an infinite sequence of semi-simple roots exists for this 
(rational) value of c. 



Figure 1: Purely imaginary roots = ik^ versus m < m for n = 1, c = -4g (left) and for n = 1, 
c = | (right). Blue dots correspond to fe+ and red dots correspond to 

Remark 2 When c = 0, the formalism of spatial dynamical systems fails. Indeed, the linearized 
system (|2.4p for e = and c = has a set of semi-simple eigenvalues of multiplicity two at k = ±ik, 
i £ Z. However, all non-zero eigenvalues are equivalent to the zero eigenvalue due to existence of 

infinitely many symmetries for the system (|2.4[) : i/j m (y) i— > r ^) m j rrn i(y)e~^ L , Vm' G Z'. Therefore, the 
behavior of all modes with m £ Z'\{n, — n} repeat the behavior of the two resonant modes with m = in 
and only these two modes are relevant for existence of localized solutions for e / and c = 0. The 
failure of the spatial dynamics formulation is related to the fact that y = x for c = 0, such that the 
second-order ODE (jl.2p can not be replaced by the PDE (|2.2p and hence it can not be written as a 
system of infinitely many second-order ODEs (12. 4p . 

Lemma 3 Consider a linear inhomogeneous equation 

d 2 + i(m - c)d y H J ip m (y) = -F m (y), Vm > m , (2.9) 
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where F m G C°(R). There exists a unique solution ip m G C°(]R), suc/i i/iai ||V ; m|lc°(R) — C|l-^mllc°(R) 
for some C > 0. 

Proof. Let VVi = e~ iay Pmiy) with a = and rewrite the linear equation (|2.9p in the equivalent 
form 

(/3 2 - 5, 2 ) ^fo) = F m (y)e^, /? 2 = 2cm ~f - - > 0. 

Since solutions of the homogeneous equation are exponentially decaying and growing as ip m ~ e^ y , 
there exists a unique bounded solution of the inhomogeneous equation in the integral form 

Vm{y) = r-= / e-P\y~y \F m {y')e ia y dy', 
such that ||<p m ||cO( R ) < ^||^m||cO(K)- D 

Remark 3 By using Lemma [3] and the Implicit Function Theorem in suitable vector spaces, one can 
solve equations of the system (I2.4p for m > mo and parameterize all components ip m (y) with m > mo 
by bounded components ip m {y) with m < mo for sufficiently small e. However, we do not perform this 
elimination at this stage, since we are going to rewrite the system (12, 4p as a Hamiltonian dynamical 
system and use a formalism of near- identity transformations and normal forms, which works easier if 
the symplectic structure of the Hamiltonian system is local. 



3 Hamiltonian formalism and normal coordinates 



We rewrite the system of second-order equations fj2.4|) as the system of first-order equations which 

2 

admits a symplectic Hamiltonian structure. Let u> = "\ + e^, where n G N and O is a free parameter. 
Let 4> m (y) = i'miy) ~ §( c ~~ m )V'm(2/) for all m G ll . The system (|2.4p is equivalent to the first-order 
system 



dy 

dtpm 

dy 



4>m + |(c - m)lp m 

-\ (n 2 + c 2 - 2cm) ip m + |(c - m) 



' f m— m\Wm\ 



(3.1) 



m—m\—m2 • 



Let bolded symbol ■?/> denote a vector consisting of elements of the set {V'mjmeZ'- The variables 
{?/>, 0, "0, 0} are canonical and the system (|3.ip is equivalent to the Hamilton's equations of motion 



dipt, 



dH 



dH 



dy d(f) m dy dtpm 
where H = H(if), cf),ip,(f)) is the Hamiltonian function given by 



m G 



(3.2) 



H 



E 



U 2 + \{n 2 + < :2 



2cm)\i> m f + -(c - m)(t/j m (j) m - ii m <t>m) + e^lVv, 



T 



E E E ^mx^-m^ 
m&L 1 mieZ 1 m 2 & 



m—m\ —rri2 



(3.3) 



S 



Let us define the discrete weighted Z 2 -space by its norm 

Vu G Z 2 (Z') : ||u||| = ^ (1 + m 2 ) s \u m \ 2 < oo. (3.4) 

meZ' 

Since ^(Z') is Banach algebra for s > ^, the convolution sums in the nonlinear system (|3.ip are closed 
if W G / 2 (Z) (Assumption [Q) and ?/> G £ 2 (Z') for s > ^. Due to the unbounded linear part, the vector 
field of the system (|3.1|) map a domain in D and to a range in X, where D and X are given by 

D = {(0, 0, 0, 0) G / 2 +1 (Z', C 4 )} , X = {(0, 0, 0, 0) G Z 2 (Z', C 4 )} , (3.5) 

for any s > |- We note that D C X and that X can be chosen as the phase space of the Hamiltonian 
system (j3.2[> . 

If W(— = W(a;), Vx G R (Assumption [T]), then W2m = w-2m, Vm G Z. In this case, the 
Hamiltonian system (|3.ip is reversible and its solutions are invariant under the transformation 

*l>(v) V'(-y)) ^ (3-6) 

In addition, the Hamiltonian function (|3.3|) is invariant with respect to the gauge transformation 

^( y ) ^ e ia ^(y), <t>{y)^e ia ct>{y), Va G R. (3.7) 

Let us define the domain for reversible solutions by 

D r = {(i>,<t>^,4>)eD: rl>{-y)=ij>{y\ 4>{-y) = -0(y)} . (3.8) 

If a local solution (ij)(y), 4>{y)) °f the system (13.1|) is constructed on y G R+ and it intersects at y = 
with the reversibility constraint 

E r = {(V>, 0,^,0) G Z? : Iim/> = 0, Re0 = 0} , (3.9) 

then the solution (t/)(y), <f>(y)) is extended to a global reversible solution in D r on y G R by using the 
reversibility transformation (I3.6H . The global reversible solution does not have an arbitrary parameter 
a induced by the gauge transformation (|3.7p . 

In order to construct non-trivial bounded solutions of the Hamiltonian system (|3.ip . we shall in- 
troduce normal coordinates for the infinite-dimensional stable, unstable and center manifolds of the 
linearized Hamiltonian system (|3.ip at the zero solution for e = (Lemma [1]). Let U_ = {m G Z' : 
m < m o} an d Z^ = {m G Z' : m > mo}. For the center manifold of the linearized Hamiltonian 
system, we set 

mGZ'_: Vm= 4 ^ 2 (y)+ /^ L, ^ = ^V / n 2 + c2-2cm[c+(y)-c-(y)]. (3.10) 
\m + cr — zcm ^ 

For the stable and unstable manifolds of the linearized Hamiltonian system, we set 

mGZ' + : i; m = ^ + f M=, 4> m = \^2cm - n 2 - c 2 [c+ (y) - c" (y)] . (3.11) 
v 2cm — n z — c z ^ 
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By using the normal coordinates (|3.10p - (|3.11|) . the Hamiltonian function (|3.3|) transforms to the new 
form 

H ^ (^ml^ml k m \c m \ ) + ^ {^m^m^m ~ ^m^m^rn) (3.12) 

m&Z'_ m£Z' + 

^~ ^ ] ^ 1 ^ ] 5171,1111,1)12(^1 ~l~ C mi)(C— ?T72 ^—7772 )( C 7Tl — 7Til— 7772 C ?77 — 777 1 —7772 ) i C m ^777 ) J 

meZ' mi ez' 7772 ez 



where 



_i_ c — m± \/n 2 + c 2 — 2cm . 
k m = , mGZ_, 



± i[c — m) ± \[2cm — — c , 
K m = - , m G Z + , 



2 _ „2 



"^777 — 7771 ™/ 

w m , mi = — 7== = = = ,- , m,mi€Z 



y^re 2 + c 2 — 2cm| |n 2 + c 2 — 2ctoi | ' 



and 

5777,7771 ,7772 



\/|n 2 + c 2 — 2cmi||n 2 + c 2 + 2cm2||ra 2 + c 2 — 2c(m — m\ — m^) ||n 2 + c 2 — 2cm\ 



for all m, mi, 1712 G Z'. The quadratic part of the Hamiltonian function in (|3.12|) for e = is diagonal 
in normal coordinates for m G Z'_ and it is block-diagonal for m E Z^_ . The Hamiltonian equations of 
motions (13.21) transform in new canonical variables to the form 



dr + 


.dH 


dr~ 

777 _ 


.dH 


dy 




dy 


dc^n 


dc+ 


dH 




dH 


dy 




dy 


dCm' 



m E Z'_ (3.13) 

ay OCm wy OCm 

and 

I 4- • 1 7 7 7 — -1 7 7 

m G Z' . (3.14) 



Since the convolution sums on {c^} me ^' involve decaying weights as \m\ — > oo and the linear unbounded 
part on {c m } me x' is linear in |m| (&;„ and grow linearly as \m\ — > oo), the vector field of the system 
(|3.13p - (|3,14p has a modified domain D' and range X' , which are given by 

£>' = {(c+ J c-,c+ J c-) € # +1 (Z',C 4 )} , X' = {(c+ c" c+ C") G #(Z', C 4 )} , (3.15) 

for s' = s — 4. (Here we have used the fact that the convolution sum acts on {ipm}m&, where ip m 
is given by in the transformation (|3.10|) - (|3.11|) .) If the space X is Banach algebra for s > |, the 
space X' is Banach algebra with respect to the decaying weights for s' > 4. The domain for reversible 
solutions become now 

^^(c+.c-.c+.c-)^': t£(-y) = c£(y), m G Z'_, c± (-y) = c^(y), m E Z' + , } (3.16) 

and the reversibility constraint becomes 

S' r = {(c + , c~, c + , c~) G D' : Imc^j = 0, m G Z'_, Rec^ = Rec~, Imc+ = — Imc~, m G Z' + } . 

(3-17) 
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In what follows, we are concerned with the reversible homoclinic orbit of the Hamilton's equations of 
motion (|3.13p - (|3.14p in normal coordinates (|3.10p - (|3.1ip . We will suppress the explicit dependence of 
H from (c + ,c~) for briefness of notations. 

Remark 4 Since > for sufficiently large negative m G the center manifold E£ © E~ is 
spanned by an infinite set of modes with positive and negative energies. Therefore, we can not use the 
technique of [5] which relies on the fact that the quadratic part of the Hamiltonian function is positive- 
definite for the non-bifurcating modes of center manifold. We can however use the technique of [6] 
which relies on the separation of slow motion for the modes, which correspond to the zero eigenvalue 
(modes and cZ n ), and the fast motion of the other modes, which correspond to non-zero purely 
imaginary or complex eigenvalues (modes for all other m). Moreover, we can simplify the technique 
of [B] by incorporating the Hamiltonian structure (|3.13|) - (|3.14p with the Hamiltonian function (|3.12p . 



4 Normal form and persistence of homoclinic orbits 

We first show that the formal truncation of the Hamiltonian function (|3.12p at the two bifurcating modes 
c+ and cZ n , which correspond to the double zero eigenvalue of the linearized system at the zero solution 
for e = 0, leads to the coupled-mode system (|1.5p for e ^ 0. Then, we derive an extended coupled-mode 
system for the modes c+ and cZ n by using near-identity transformations of the Hamiltonian function 
and prove persistence of a reversible homoclinic orbit in the extended coupled-mode system. We will 
assume from now on that < c < n. 

Let us consider the subspace of the phase space of the Hamiltonian system (I3.13p - (|3.14[) : 

S={ c +=0,VmeZ'\{n}, c~ = 0, Vm G Z'\{-n}} . (4.1) 

Denote Q = S^. If the Hamiltonian function H is formally constrained on the subspace S, the 
expression (|3.12p takes the form 



H\ 



^I4| 2 , ^l c -n! 2 _ W 2n (c+C_ n + C+C_J 

n - c n + c yjn 2 - c 2 




Since (n, —re) G , we use the symplectic structure (|3.13p to generate a system of first-order ODEs for 
normal coordinates c+ and cZ n - By using the new independent variable Y = ey and the new dependent 

+ c~ 

variables a = ^— and b = —^==. we recover the ODE system 



i(n — c)a' + £la = W2 n b + cr(\a\ 2 + 2\b\ 2 )a, 
-i(n + c)b' + nb = w 2n a + a(2\a\ 2 + \b\ 2 )b, 



2 , (4-2) 



where the derivatives are taken with respect to Y = ey. The system fj4.2|) is nothing but the coupled- 
mode system f)1.5|) after the separation of the traveling variable Y = X — cT and the stationary variable 
T in the transformation 



a{X, T) ^ a(Y)e- lQT , b(X, T) ^ b(Y)e~ inT ', Y = X - cT, 

as 

(a homoclinic orbit) (see Section 1) for W2 n an d |^| < ^o> where = |^2n| — 



with the correspondence Y = \/n 2 — c 2 £ and f2 = C J . The system (|4,2p has a localized solution 
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The ODE system (|4,2|) is invariant with respect to translation a(Y) — > a(Y"— Yo), i>(Y") — ► 6(Y"— lo) for 
all Yo £M and gauge transformation a(Y") — > e ia a(Y~— Yo), 6(Y") — > e ia 6(Y"— Yo) for all oel. Therefore, 
any solution of the system is continued with a two-parameter group of symmetry transformations. 
However, these parameters are set uniquely in the reversible homoclinic orbit of Definition [TJ In 
addition, we note that although the ODE system (|4.2p is formulated in the four-dimensional phase 
space, it has two conserved quantities on Y E R related to the translational and gauge symmetries. 
Indeed, the Hamiltonian H\$ and the quadratic function 

Q = |c+| 2 - \cZ n \ 2 = {n- c)\a\ 2 - (n + c)|6| 2 (4.3) 

are constants on Y E E. As a result, localized solutions of the system (|4.2|) are defined on a subspace 
which obeys a planar Hamiltonian system. This planar Hamiltonian system is given by the second 
equation of the system (|1.1U|) in variables (cp,4>). We recall that reversible homoclinic orbits of planar 
Hamiltonian systems are structurally stable with respect to parameter continuations. 

To incorporate the ideas of integrability of the Hamiltonian coupled-mode system on the subspace 
S and persistence of the reversible homoclinic orbit in the planar Hamiltonian system, we extend the 
coupled-mode system by using near-identity transformations and the normal form theory. 

Lemma 4 Let < c < n, such that n2 + c2 1/ . p or each N G N and sufficiently small e, there is 
a near-identity, analytic, symplectic change of coordinates in a neighborhood of the origin in X' in 
13.15\) . such that the Hamiltonian function H in h3.12\i transforms to the normal form up to the order 
ofO(e N+1 ), 



H ^ 1 (km\ c m\ k m \c m \ ) + ^ ^ iyK m C m C rn K m C m C m ^ 



+ tH s (4, cZ n ) + eH Q (c+,cZ n , c+ c") + e N+1 H R (c+, cZ n , c+, c"), (4.4) 

where Hs is a polynomial of the degree 2N + 2 in (c^,cZ n ) on S, Hq is a polynomial of the degree 2N 
in (c^,cZ n ) on S and of the degree 4 in (c + ,c~) on Q with no linear terms with respect to (c + ,c~) on 
Q, and Hr is a polynomial of the degree 8N + 4 in (c+, cZ n ) on S and of the degree 4 in ( c+ 5 c~) on 
Q. All components Hs, Hq and Hr depend on e, such that H$ and Hq are polynomials in e of the 
degree N — 1 and Hr is a polynomial in e of the degree 3N — 1. The reversibility §3. 6\) and gauge {3. 7\ ) 
transformations are preserved by the change of the variables. 

Proof. The existence of a near-identity symplectic transformation that maps H to the form (|4.4p 
follows from the fact that the non-resonance conditions Ikq — are satisfied for any / G Z and 

m G U since k,q = and all eigenvalues are semi-simple for e = 0. See [5] for an iterative sequence of 
symplectic transformations. The transformation is analytic in a local neighborhood of the origin in X' 
as the vector field of the Hamiltonian system f)3.13p — (|3.14j) is analytic (given by a cubic polynomial). 
The reversibility f|3.6|) and gauge f|3.7[) symmetries are preserved by the symplectic change of variables 
[5]. The count of the degree of polynomials Hs, Hq and Hr follows from the fact that the vector field of 
the Hamiltonian system (|3.13p - (|3.14p contains only linear and cubic terms in normal coordinates, while 
the near-identity transformation of (c + ,c~) on Q up to the order of 0(e N+1 ) involves a polynomial in 
e of the degree ./V and a polynomial in (c+, cZ n ) on S of the degree 2iV + 1. □ 

Example 2 For N = 1, the near-identity transformation for c^, m G Z'_\{ra} takes the explicit form 

— < ~"m T^p [^m,?iC n + W m ^— n C_ n + G (^gn,n,n{c n ) C_ n 5 m ^ n + g— n ^— n — r \c n \ C_ n 5 m ._ ra ) j + 0(6 ), 
Km 
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where c+ is a new variable. Similar explicit formulas can be obtained for c m , m £ n} and for 

c*, m € It' + . We note that the tilde signs are omitted for new variables in the expression (|4.4p . 

Remark 5 For each N € N, the subspace S 1 defined by ()4. 1[) is an invariant subspace of the Hamiltonian 
system f|3.13j) — (j3.14[) with the Hamiltonian function (|4.4p truncated at -ff^j = 0. The dynamics on S" is 
given by the four-dimensional Hamiltonian system 

d4 _ dHs dc^ _ dH^ 

dY ~ dct 1 dY ~ 'dctj [ } 

where Y = ey. If N = 1, the system (I4.5P transforms to the coupled-mode system (14.21) in variables 



a = — £s= and b = -, " . If JV > 1, this system is referred to as the extended coup led- mode system. 

yn- c v n + c 

Lemma 5 Lei 7^ for a given n £ N. For eac/i IV £ N and sufficiently small e, there exists a 

1 — ^ 2~ 

reversible homoclinic orbit of the system ft4-5\ ) for \Q\ < = \ w2n\ ra ~ c • Moreover, the solution for 
the homoclinic orbit satisfies the global bound 

\c+(y)\ < C+e-^l, |C n (y)| < CLe"^, Vy G R, (4.6) 

/or some e-independent constants 7 > and C-j- > 0. 

Proof. Due to the gauge- invariance of the polynomial Hamiltonian function H$ = Hs(c^, c~ n ), it must 
satisfy the partial differential equation [2]: 



±Hs(e ia c+,e ia cZ 



_ + 8H S _+dH s dH s 8H S _ 

C n 1 C n — L T C_ n _ C_ n _ U. ' ) 

a=0 OCn OCj OC_ n OC_ n 



It follows from the system (|4.5p and the relation (|4.7p that Q = |c+| 2 — |c_ ra | 2 is constant in Y € M. [2\. 
If localized solutions exist, then Q = 0. Let us represent the solutions in the general form 

c+ = Vp + qe^ +w , cZ n = Vp^qe^- id , (4.8) 

where (p,q,9,(p) are new real- valued variables. Using the chain rule for Hs = Hs(p,q,9,ip), we find 
that (p is a cyclic variable of the Hamiltonian function Hs and q is constant due to the gauge invariance 
(|4.7p . Setting g = 0, we find that (/?, 0) satisfy a planar Hamiltonian system, while <p is found from a 
linear inhomogeneous equation: 

n dp _ 8H S] 9 ^_^£| 9^_^i f 4 q\ 

dY ~ 89 ' 9= °' dY ~ dp ' 9=0 ' dY " 8q l9=0 ' 1 j 

where Hs\ q =o is a function of (p,9). If e = 0, the system (|4.9p reduces to the ODE system (|1.10p 
rewritten in new coordinates and for cj) = ^fpe % ® . The vector field of the extended coupled-mode system 
(|4.5|) is given by polynomials in c+ and c~ n of the degree 2iV-|-l and in e of the degree N—l. Recall that 
the coupled-mode system fj4.2|) admits a reversible homoclinic orbit for W2 n 7^ and \Q\ < Qq, where 
Slo = \u)2n\ ^^r" ~- Since a reversible homoclinic orbit is structurally stable in a planar Hamiltonian 
system with an analytic vector field, the homoclinic orbit persists in the extended coupled-mode system 
(|4.5p for sufficiently small e. □ 
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5 Construction of local center and center-stable manifolds 



We study solutions of the Hamiltonian system of equations (|3.13p - (|3.14p after the normal- form transfor- 
mation of LemmaHl We construct a local solution on y E [0, L/e N+1 ] for some e-independent constant 
L > 0, which is close to the homoclinic orbit of Lemma by the distance Ce N for some e-independent 
constant C > in vector norm on X' . This solution represents an infinite-dimensional local center- 
stable manifold and it is spanned by the small bounded oscillatory and small exponentially decaying 
solutions near the exponentially decaying homoclinic solution with the decay bound (|4,6p . Parameters 
of the local center-stable manifold are chosen to ensure that the manifold intersects at y = with the 
symmetric section . This construction completes the proof of the main Theorem [2j 

By using Lemma [J] and the explicit representation of the Hamiltonian function (|4.4p , we rewrite the 
Hamiltonian system of equations in the separated form 



dy 

dc -n 

dy 
dc 
dy 



/i+ c + + ^ + c :„, + F+( c +, c :„) + F+(4,c: n ,c)j + e JV+i F+( c +, c ;, c ), (5.1) 

(j,Jc+ + v-cZ n + Fg (c+,cZ n )+Fn(c+,cZ n ,c) + € N+1 F R (c+,c~,c), (5.2) 



A e c + eF Q (c+, c_ n , c) + e N+1 F R (c+,c n , c), 



(5.3) 



where c denotes all components of (c + , c , c + , c~) in Q = S , (/i^,^,A e ) denote the coefficient 
matrix for the linear part of the system and (F s , Fq , F R , Fq, Fr) denote the nonlinear (polynomial) 
part of the system. For briefness of notations, we do not rewrite the subsystem (|5.ip ~ (|5.2p for variables 
(c+, cl n ) and we do not write dependence of the nonlinear functions from these variables. The variables 
c are equivalent to tilde-variables in Example [2] after the near-identity transformations, but the tilde- 
notations are dropped for simplicity of notations. 



Lemma 6 Let W(x) satisfy Assumption{]\ and W2 n 7^ for a given n £ N. For sufficiently small e, 
the linearized system 115. 1\) ~ [573\) at the zero solution for \Q\ < fi = \ w 2n\ w ~ c an d < c < n is 
topologically equivalent to the one for e = 0, except that the double zero eigenvalue of the subsystem 
115. ij) -( f5Jg]) splits into a pair of complex eigenvalues to the left and right half-planes. 

Proof. Since all non-zero eigenvalues of the linearized Hamiltonian system (I5.3P at the zero solution 
are semi-simple at e = 0, they are structurally stable in the perturbation theory for sufficiently small 
e 7^ 0. The matrix operator A e is a polynomial in e and Ao is a diagonal unbounded matrix operator 
which consists of (ik^, —ik^, —ik^ik^) for m £ Z'_ and of (k+ , k~ , , k~ ) for m G Z' + . The matrix 
operator with elements ew7 m ,m' represents a small perturbation to Ao if the vector of Fourier coefficients 
W is in Z X (Z) C Z 2 (Z) for s > ^ according to Assumption [TJ 

The coefficients f£~ and vf are polynomials in e and 

_l_ in _ iW2n + —iw2n - 

At e = 0, the linearized subsystem (|5.1|) - ([5.2|) corresponds to the linearized coupled-mode system (|4.2p . 
Its characteristic equation is given by 

(n 2 - c 2 )k 2 - 2iecttK + e 2 (fl 2 - w% n ) = 0, 
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with two roots 



K± 



iQc ± \J (re 2 



2n 



n 



2 ^ 2 



- c z 

\w2r, 



The two roots have ReK+ ^ if \Q\ < Oq, where fio = \w2n \ s ^ a "~ c " ■ Under the same assumption on 
W, perturbation terms in fif and vf- are small compared to the leading-order terms and f^, such 
that the pair persists in the left half-plane and right half-plane of the K-plane. □ 



Corollary 1 For sufficiently small e, a local neighborhood of the zero point in the phase space X' can be 
decomposed into the subspaces determined by the spectrum of the linearized system at the zero solution 

x' = x h ex c ex u ex s , (5.4) 

where X s and X u are associated to the subspaces E s and E u of Lemma [71 while X n and X c are 
associated to the subspaces E c+ © E c on S and Q respectively. 

By Remark [5l the truncated system (|5.ip - (|5.3p with F^ = and F# = admits an invariant 
reduction on S. By Lemma (H the extended coupled-mode system (14. 5|) on S has a reversible homoclinic 
orbit which satisfies the decay bound (|4.6p . This construction enables us to represent the solution of 
the subsystem (I5.1|) - (15.2p and its complex conjugate in the form [c+ , cZ n , c+ , cZ n ] T = co(ey) + c n (y), 
where co(ey) is the homoclinic orbit of Lemma [5] and c n (y) is a perturbation term. By using the 
decomposition, we rewrite the system (|5.ip - (|5.3p in the equivalent form 



dc h 
dy 
r/c 
dy 



eh h {c G )c h + eG Q (c )(c h , c) + e N+1 G R (c + c h , c) 
A e c + eF Q (c + c h , c) + e N+1 F R (c + c h , c), 



(5.5) 
(5.6) 



where A^(co) is a 4-by-4 linearization matrix of the extended coupled-mode system (14. 5p and its complex 
conjugate at the solution co(ey) and (Gq,Gr) denote the nonlinear part of the subsystem (|5.ip - (|5.2p 



and its conjugate. We note that the function Gq combines nonlinear terms in c n from the functions F s 
and the nonlinear terms in c from the functions Fq . For simplicity of notations, we say that c n G X) 
and c G Xj- in the decomposition X' = X^ © Xj- . 



Example 3 For N = 1, the linearization matrix A/j(cq) takes the explicit form 



n—c 



(n—c) 2 n 2 — c 2 
2n 2 f 75+cl n 
t=? i „2 „2 



"'2n 



2<jc„ 



n+ 2 



\J1. 



n+c 



-+2 
{n—c) 2 
2actc~ 



, 2u\c+\ 2 

~T~ — 1 T~ 

2aci c~ 



+ 



(n+c) 2 



i 2 — i 



n—c 



(n-c) 2 

2actcZ n 
— 5 ? — 

Met? 



+ (n-c) 2 + 



._ 



2<rctc_ n 
5 3 — 

n 4 — e 4 

(n+c) 2 

2<tc^~c~ 

+ 



(n+c) 2 



n z — c^ 



Vn 2 -c 2 ' « 2 -c 2 
Q _ 20-lc+l 2 _ 2 CT |c:J 2 
n+c n 2 — c 2 (n+c) 2 



where components of [c+, c_ n , c^, c_ n ] T are evaluated at the solution Co of the extended coupled-mode 
equation (|4.5p . As \y\ — > oo, the matrix A/ 1 (cq) converges to the form 



A„(0) 





««2n 


n—c 


\/n 2 — c 2 


»«2n 


Q_ 


n 2 — c 2 


n+c 































n—c 


\fn 2 — c 2 




n 


n 2 — c 2 


n+c 
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such that || A/j,(co) — Af l (0)\\ Xhl ^, Xh < Ce n ^ for some C > and 7 > according to the decay bound 

ESI). 



/ 2 2 

Lemma 7 Lei W2n 7^ for a given n G N, |Q| < = l w 2n| n ~ c > and < c < n. Consider the 
linear inhomogeneous equation 

^-eA h (c )c h = F h (y), (5.7) 
dy 

where G C°(R). T/ie homogeneous equation has a two-dimensional stable manifold spanned by the 
two fundamental solutions 

s i = c o(2/)> s 2 = o-ico(y), (5.8) 
where o~\ is a diagonal matrix of (1, 1, —1, —1). // components o/F/ l (y) satisfies the constraints Vy G R 

(F ft ) 1 ( y ) = (F h )i(-y), (F ft ) 2 (y) = (F&) 2 (-y), (F fc ) 3 (y) = -(*\)i(y), (F h ) 4 (y) = -(F /l ) 4 (y), 

(5-9) 

i/ien i/iere exists a two-parameter family of solutions Ch G C?(H.) in i/ie /orm = aisi(y) + «2S2(y) + 
Ch(y), where (01,02) o- r ^ parameters and Ch(y) is a particular solution of the inhomogeneous equation 
|5, 7| ), such that ||c7,,||c7°(]R) < y[|F&||c°(tr) f or an independent constant C. 

Proof. The existence of the two-dimensional kernel (|5.8|) follows from symmetries of the extended 
coupled-mode system (|4.5p with respect to translation and gauge transformation. Since the subspace 
Xh associated with A^(0) is invariant under A/j(cq), the kernel is exactly two-dimensional and the 
other two fundamental solutions of the homogeneous equation are exponentially growing. The adjoint 
homogeneous equation ^ = — eA^(co)u has also a two-dimensional stable manifold spanned by the 
two fundamental solutions 

s* 1 = a 2 c' (y), s* 2 = fJ 3 co(y), (5.10) 

where 02 and 03 are diagonal matrices of (1, — 1,-1,1) and (1,-1,1,-1) respectively. Unless the vector 
function F/j is orthogonal to {s*^}, a solution of the linear inhomogeneous equation (|5.7p grows 
exponentially as |y| — > 00. However, (s*,F^,) = (si^F/J = if the constraints (|5.9p are satisfied. By 
the Fredholm theory, F^ is in the range of the linear unbounded operator ^ — eA/ l (co), such that there 
exists a solution G C®(M) of the inhomogeneous equation (|5.7p such that Hc/Jl^o^ < 7-||F^|| C 0( K ) 
for an e-independent constant C. A general solution of the inhomogeneous problem has the form 
Ch = a±si(y) + «2S2(y) + Ch{y)> where (a±, 02) are parameters. □ 

Lemma 8 The nonlinear part of the vector field of the system \5. 5\) - [57S\) is bounded in a local neigh- 
borhood of the zero point in X' = X^ © X^ by 



||Gfl||x fc < N R (J|c + c h \\ Xh + \\c\\ x jl) , ||F fi || x ± < M R (J|c + c h \\ Xh + ||c|| x xj , (5.11) 
\\G Q \\ Xh < N Q (\\c h \\ 2 Xh + ||c||^) , \\Fq\\ x > < M Q (||c + c h \\ Xh + \\c\\ x ±^ \\c\\ x ±, (5.12) 
for some N R , Mr, Nq, Mq > 0. 



Proof. The system (|5.5p - (|5.6p is semi-linear with polynomial vector field for a finite N G N defined on 
the domain D' C X' , where D' and X' are given in (|3.15p . If X is the Banach algebra for s > 5, then 
X' is Banach algebra for s' > j. The derivation of estimates (|5,lip - (|5.12p follows similarly to Lemmas 
2 and 3 in [5] ■ The characterization of Gq and Fq is based on the fact that the Hamiltonian function 
Hq is quadratic with respect to c by Lemma □ 
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Remark 6 Using Corollary [Q we denote A c = A e |x c , A u = A e |x u , A s = A e |x s for a block-diagonal 
decomposition of A e on the invariant subspaces X c , X u and X s respectively. We also denote the coor- 
dinates of the decomposition by c = [c c ,c n ,c s ] and the projection operators by P c , P u , P s respectively. 

Theorem 3 (Local center-stable manifold) Let a G X c , b G X s and (ai,a 2 ) £ C 2 be small such 
that 

\Hxc < C a e N , \\b\\ Xs < C b e N , \ ai \ + \a 2 \ < C a e N . (5.13) 

for some e-independent constants C a ,Cb,C a > 0. Under the conditions of Lemma® there exists 
a family Ch(y; a, b, a\, a 2 ) and c(y; a, b, a±, a 2 ) of local solutions of the system Ii5.5\ )- [5l)\) such that 
c c (0) = a, c s = e yAa b + c s (y), c h = aiSi(y) + a 2 s 2 (y) + c h (y) with uniquely defined c s (y) and c h (y), 
and the local solutions satisfy the bound 

sup \\c h (y)\\ Xh < C h e N , sup \\c(y)\\ x x < Ce N , (5.14) 

y£[Q,L/e N + 1 ] ye[0,L/e N + 1 ] h 

for some e-independent constants L > and Ch,C > 0. 



Proof. We modify the system (|5.5p - (|5.6p by the following trick. We multiply the nonlinear vector field 
of the subsystem (]5.6f) by the cut-off function X[a,y ]{v)i such that 

dc 

— = A e c + ex[o, yo ](y)F Q (co + c h , c) + e N+1 X[o, yo ](y)F R (c + c h ,c), (5.15) 

where X[o, yo ](y) = 1 for y G [0, yo] C M for some yo > and X[o,y ]( y ) = otherwise. Similarly, we 
multiply the nonlinear vector field of the subsystem (|5.5p by the cut-off function X[o,yo]( v ) an< ^ a ^ 
symmetrically reflected vector field multiplied by the cut-off function X[- yo ,o](y)i such that 

^-eA h (c )c h = ex[o, yo ](y)G Q (c )(c h ,c) +e N+1 X[o, yo ](y)G R (c + c h ,c) 

£X[-y ,o](y)G* Q (c )(c h ,c) + e N+1 X[^y ,o](y)G* R (c + c h ,c), (5.16) 

where (G* QR ) lj2 (y) = (G QtR ) lj2 (-y) and ^Q R hM = -(Gq^) 1,2(3/) on y G [-y ,0). We are 
looking for a global solution of the system (|5.15|) - (|5.16|) in the space of bounded continuous functions 
C£(R). This global solution on y G M corresponds to a local solution of the system (|5.5p - (|5.6p on the 
interval y G [0, yo] C R. 

The imaginary axis lies in the resolvent set of A u and A s and 

for some -Ko > 0. Let c s (y) = e^b + c s (y) and look for solution c s (y) and c u (y) of the system 
(|5.15p projected to X s and X u with operators P s and P u . By Lemmas El El and the Implicit Function 
Theorem, there exists a unique map from C°(1R, X^ © X c ) to C°(R, X M © X s ) parameterized by b such 
that 

sup [||c u (y)|| Xa + ||c s (y)||x 3 ] < ||b|| x .+eMi sup [(1 + \\c h {y)\\ Xh + ||c c (y)|UJ ||c c (y)|UJ 

y6[0,J/o] l/6[0,s/o] 

+e N+l M 2 sup [l + ||c fc (y)||x h + IMl/)lkJ , (5-18) 
2/e[0,y o ] 
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for some M 1 ,M 2 > 0. 

Let Ch(y) = «iSi(y) + a 2 s 2 (y) + Ch(y) and look for solution Cft(y) of the system (|5.16p . By the 
Hamiltonian structure of the system (|5.5p . the vector field satisfies the constraints (GQ^^^y) = 
— (GcQ t ji)i t 2(y) on y G [0, yo]- By the construction of the modified vector field, it satisfies the constraint 
(|5.9p . By Lemmas [3 [8j the bound (|5.18p . and the Implicit Function Theorem, there exists a unique 
map from CP(R, X c ) to CP(R, X7J parameterized by (ai,a2) and b such that 

sup \\c h {y)\\ Xh < \ai\ + \a 2 \ + M 3 sup ||c c (y)||^ c + e JV M 4 sup [1 + ||c c (y)|| x J , (5.19) 

!/e[0,j/o] ?/6[0,2/o] 2/6[0,yo] 

for some M 3 ,M 4 > 0. 

Since the spectrum of A c consists of pairs of semi-simple purely imaginary eigenvalues, the operator A c 
generates a strongly continuous group e y!Vc for any y 6 R on X c such that 



sup || e 



wAe lk,~x c <tf, (5.20) 

for some -ftT > 0. By variation of constant formula, the solution of the system ()5. 15j) projected to X c 
can be rewritten in the integral form 

c c (y) = e^a + e [ V e(y~y'^P c [F Q (c (ey') + c h (y'), c(y')) + e N F R (c (ey') + c h (y'), c(y'))] dy' , 
■/ o 

(5.21) 

where a = c c (0). By using the bound ()4.6p for co(ey) and the bounds (15.18P and (|5.19p on the 
components c njS and c^, we derive from the integral equation (|5.2ip that 



/ ryo 
sup \\c c (y)\\ Xc < K [\\&\\x c + \h\\x s + \ai\ + \a 2 \+eM 5 / \\c (y)\\ Xh \\c c (y)\\x c dy 
e[0,yo] V Jo 

ryo \ 
+ey M 6 sup ||c c (y)||| c + e N+1 M 7 / ||c (y)||x^y + ^l/oMg sup ||c c (y)|| x l 



y£[o,j/o] Jo ye[o,j/o] 

for some M§, Mq, M-?, M$ > 0. By the Gronwall's inequality, we have thus obtained that 

sup ||c c (y)|| Xc < Ke eKM ^o° \My)\\x h d y + + ^ + ^ 

ye[o,j/ ] 

+e JV M 9 + ey M 6 sup ||c c (y)||^ c + e^+^Mg sup ||c c (y)||* c ) , (5.22) 
y6[0,y ] ye[0,yo] J 

for some Mg > 0. Here we can use the decay bound (|4.6|) which implies that e \\co(y)\\x h dy < C 
for some e-independent C > 0. By using the same bound for the exponent and letting yo = L/e M , we 
can see that we can choose M < iV + 1, where the value M = N + 1 gives the balance of all terms in 
the upper bound (|5,22p . If arbitrary vectors a, b and (a\,a 2 ) satisfies the bound (|5,13p . then we have 
constructed a local solution c c (y) which satisfies the bound 

sup \\c c (y)\\ Xc < C c e N (5.23) 

2/S[0,L/e JV + 1 ] 

for some C c > 0. By using the bounds (I5.13p . (|5.18j) . (|5.19j) . and (I5.23p . we have proved the bound 
(|5.14p for some e-independent constants Ch, C > 0. □ 
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Remark 7 One can prove Theorem by using the contraction mapping principle and the integral 
formulation for the local center manifold of the system (|5.5p - (|5.6p . This approach was undertaken in 
Section 4 of [5] (see their Theorem 4). We have avoided this unnecessary complication with the explicit 
decomposition (|5.4p and analysis of the system (|5.5p - (|5.6p decomposed into subsystems. Similar direct 
methods of analysis have been applied to problems without Hamiltonian structures such as the quasi- 
linear wave equation in [6], where an iteration scheme was employed to prove the bound on small local 
solutions along the local center-stable manifold. 

Proof of Theorem By Theorem [31 we have constructed an infinite-dimensional continuous family of 
local bounded solutions of the system (I5,5p - (|5.6p on y G [0, L/e N+1 ] for some e-independent constant 
L > 0. The solutions are close to the reversible homoclinic orbit of the extended coupled-mode system 
(|4.5p in the sense of the bound (|1.14p . It remains to extend the local solution to the symmetric interval 
y G [— L/e N+l , L/e N+1 ] as the local reversible solution with the reversibility constaints (I3.16p . To do 
so, we shall consider the intersections of the local invariant manifold of the system (|5.5p - (|5.6p with the 
symmetric section T,' r defined by (|3.17H . 

Since the initial data c c (0) = a in the local center-stable manifold of Theorem [3] are arbitrary, the 
components of a can be chosen to lie in the symmetric section T,' r , such that 

Im(a)+ = 0, Vm G Im(a)~ = 0, Vm G Z'_\{-n}. (5.24) 

This construction still leaves infinitely many arbitrary parameters for 

Re(a)+, Vm G Z'_\{n}, Re(a)~, Vm G Z'_\{-n} (5.25) 

to be chosen in the bound (|5.13p . The initial data c^(0) and c su (0) are not arbitrary since we have 
used the Implicit Function Theorem for the mappings (|5.18p and (|5.19p . Therefore, we have to show 
that the components of b and (01,02) can be chosen uniquely so that the local center-stable manifold 
intersects at y = with the symmetric section E^.. 

We note that there are as many arbitrary parameters s and (a\, 02) in the local center-stable manifold 
as there are remaining constraints in the set Y/ r . First, let us consider constraints in the set Y/ r for 
m G namely 

Rec+(0) = Rec~ (0), Imc+(0) = -line" (0), m G Z' + . (5.26) 
Let c s = e yAs b + c s (y) and rewrite the constraints in the form 

Re6 m + Re(c s ) m (0) = Re(c u ) m (0), Im6 m + lm(c s ) m (0) = -Im(c u ) TO (0), (5.27) 

where all terms are of order O(e^) and the vectors c s (0) and c u (0) depend on b in higher orders in 
e. By the Implicit Function Theorem, there exists a unique solution of the constraints for b such that 
[|b[|x s satisfies the bound (|5.13p . 

Finally, let us consider constraints in the set Y! r for components of c^, namely 

Imc+(0) = 0, lmcl n (0) = 0. (5.28) 

Let Ch = a\S\ + «2 S 2 + c s (y) and note that si and S2 violate the constraints (|5.28p . Let a and b be 
chosen so that c(0) belongs to the set H' r . By a construction of the vector field in the system (|5.16p . if 
c(y) lies in the domain D' r of reversible solution and a% = 02 = 0, then Gq r = Gq,r on y G [— yo, yo] 
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and the global solution c s (y) constructed in Theorem intersects the set X', at y = 0. Therefore, the 
choice ai = «2 = satisfies the constraints (|5,28p identically. 

We have thus constructed a family of reversible solutions in the symmetric interval y G [— L/ e N+l , L / e N+ 
while preserving the bound (|1.14|) . Tracing the coordinate transformations used in our analysis back 
to the original variable if)(x,y), we have thus completed the proof of Theorem [2j □ 



6 Discussion 



We have proved that a moving gap soliton of the Gross-Pitaevskii equation (jl.ip with the periodic 
potential V(x) is surrounded by the oscillatory tails which are bounded on finite intervals of the spatial 
scale. Because the center manifold is infinite-dimensional with the sign-indefinite Hamiltonian function, 
we are not able to exclude the polynomial growth of the oscillatory tails in the far-field regions. This 
construction of traveling solutions on a finite spatial scale is related with the finite-time applicability 
of the coupled-mode equations 111. 5ft for the Cauchy problem associated with the Gross-Pitaevskii 
equation dTTTJ) [T2] . 

It would have been a drastic improvement to the constructed theory if we could extend the analysis 
of oscillatory tails to the infinite spatial scale by proving existence of global solutions with a single bump 
and bounded oscillatory tails. In many problems with finite-dimensional center manifolds associated 
with semi-simple purely imaginary eigenvalues, such constructions of global center-stable manifolds are 
possible and the proof of persistence of bounded solutions with oscillatory tails can be developed [3]. 

We will show that the basic evolution models for moving gap solitons in periodic potentials exhibit 
infinite-dimensional center manifolds in the spatial dynamics formulation. In particular, we can think 
of three possible generalizations of the Gross-Pitaevskii equation (jl.ip . given by the complex- valued 
Klein-Gordon equation 

Ett-E xx = V(x)E + a\E\ 2 E, (6.1) 
the regularized Gross-Pitaevskii equation 

iEt = -E xx + iE xxt + V{x)E + a\E\ 2 E, (6.2) 

and the discrete Gross-Pitaevskii equation 

iE n = —E n+ \ — E n -i + V n E n + a\E n \ 2 E n . (6.3) 

When V(x) = or V n = 0, the spectrum of the linear part of the Klein-Gordon equation (|6.ip is 
unbounded from both above and below, while that of the regularized and discrete Gross-Pitaevskii 
equations (16. 2p and (|6.3p is bounded from both above and below. We look at the traveling solutions of 
these equations in the form 

E(x, t)=J2 ^m{y)e^- tult , y = x-ct, (6.4) 

for the linear limit a = with no potential V(x) = 0. (In the case of the lattice equation (|6.3p . we 
use the traveling ansatz (|6,4p at x = n for n G N.) As a result, we obtain uncoupled linear ODEs or 
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differential advance-delay equations for amplitudes {ipm(y)}meZ' which are solved with the substitution 
V'm(y) = Vm £ Z'. All roots k are found from the following characteristic equations 



in 



2 



-icK 3 + (u — 1 + mc)K 2 + i ( c — m + mw + — k — <j H — ^-(1 — a;) = 0, (6-6) 

to — ick + 2 cosh + — ^ = 0, (6-7) 



(1 - + i(m - 2cuo)k + lo - — = 0, (6.5) 

2„\ ™2 



which correspond to the three relevant models (|6.ip - (|6.3p . It is easy to see that the characteristic equa- 
tions (|6.6p and (|6.7p have at least one purely imaginary root p for any m G Z', while the characteristic 
equation (16. 5p has two purely imaginary roots p for sufficiently large values of \m\ on m £ Z'. There- 
fore, the dimension of the center manifold associated with the linearized system at the zero solution 
for V(x) = is infinite in all three models (|6.ip - (|6.3p . 
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